X-ray Detection from Bona-fide and Candidate Brown Dwarfs in 
the p Ophiuchi Cloud with Chandra 
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ABSTRACT 

We present results of an X-ray search from bona-fide and candidate brown dwarfs in the p 
Ophiuchi cloud cores with the Chandra X-ray Observatory. The selected areas are two fields 
near the cloud center and are observed with the ACIS-I array of a 17' x 17' size and a '--^lOO ks 
exposure. Among 18 bona-fide and candidate brown dwarfs listed by the infrared spectroscopy, 
we find X-ray emission from 7 sources above 99.9% confidence level. Therefore ~40% of the 
infrared-selected brown dwarfs in this cloud emit X-rays. For the brightest 4 sources, the X-ray 
spectra are made and are fitted with a thin-thermal plasma model of a temperature 1-2.5 keV. 
The X-rays are also time variable with rapid flares from 2 of the brown dwarfs. Assuming 2 keV 
temperature and using the empirical relation of Ay vs. iVn, we estimate the X-ray luminosity 
or its upper limit of the other faint or non- X-ray sources. The X-ray luminosity {Lx) of the X- 
ray-detected sources is in the range of 0.3-90x10^^ ergs s^^, while the luminosity ratio of X-ray 
to bolometric (ijc/iboi) is lO^'^-lO^^, similar to those of low-mass pre-main-sequence and dMe 
stars. All these results suggest that the X-ray origin of brown dwarfs is the same as low-mass 
stars; strong magnetic activity at the stellar surface. 

Subject headings: stars: low-mass, brown dwarfs — X-rays: stars — ISM: individual (p Ophiuchi cloud) 



1. INTRODUCTION 

Brown dwarfs are sub-stellar objects with mass 
well below the hydrogen burning limit (~0.08 
Mq), which fill up the gap between stars and plan- 
ets. With no stable nuclear burning, the energy 
source of brown dwarfs is gravitational contrac- 
tion, hence brown dwarfs become cooler and less 
luminous with the increasing ages. ROSAT de- 
tected X-rays from several brown dwarfs and their 
candidates in star-forming regions; Chamaeleon, 
p Ophiuchi [p Oph), and Taurus (Neuhauser & 
Comreron 1998; Neuhauser et al. 1999). X- 
ray flares from field brown dwarfs were also found 
(Fleming, Giampapa, & Schmitt 2000; Rutledge 
et al. 2000). These X-ray features are simi- 
lar to those of low-mass stars, hence the X-rays 
are likely to be magnetic origin. However the 
standard (a-oj) dynamo to produce the magnetic 



activity may not present, because brown dwarfs 
are fully convective, hence there are no anchor 
points of the magnetic field in the star interior 
(Drake et al. 1996). In order to address the 
X-ray features and to study the emission mech- 
anisms of brown dwarfs, our current knowledge is 
still very poor, due mainly to the limited sensi- 
tivity of the previous instruments. This observa- 
tional constrain is largely relaxed by the Chandra 
X-ray Observatory of wide band sensitivity (0.5- 
10.0 keV) coupled with unprecedented spatial res- 
olution of ^0.5" (Weisskopf, O'dell, & van Spey- 
broeck 1996). We systematically search and study 
X-rays from brown dwarfs in one of the nearest 
star-forming regions, using the data of two deep 
ACIS-I exposures on the p Oph molecular cloud 
cores. 
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2. OBSERVATIONS AND DATA RE- 
DUCTION 

Two Chandra observations were made on the 
central region of the p Oph cloud with the ACIS- 
I array consisting of four abutted X-ray CCDs. 
The first observation (here and after, obs.l) cov- 
ered a 17.'4xl7.'4 area including cores B, C, E, 
and F, while the second observation (obs.2) cov- 
ered the center of core A (Motte, Andre, & Neri 
1998). The level 2 data are retrieved from the 
Chandra X-ray Center (CXC) archive, in which 
the data degradation caused by the increase of 
charge transfer inefficiency (CTI) in orbit is cor- 
rected. X-ray events are selected with the ASCA 
grades 0, 2, 3, 4, and 6. After the processing, 
fslOO ks effective exposure time is obtained from 
each observation. The log of the observations is 
listed in Table 1. 

In these two fields, 18 late M dwarfs have been 
reported, based on the water vapor absorption 
at A = 2.4 2.5 /im (Wilking, Greene, & Meyer 
1999; Cushing, Tokunaga, & Kobayashi 2000). 
Among them, 8 sources have the upper limit of 
mass less than 0.08 M0, which we call "bona- 
fide brown dwarfs" (here and after, BDs). The 
other 10 sources have a mass in the transition re- 
gion of 0.08 M0, hence called "candidate brown 
dwarfs" (CBDs). Their names and spectral types 
are shown in Table 2. 

3. ANALYSIS AND RESULTS 
3.1. Source Detection 

Using the wavdetect^ command, we pick up 
~100 X-ray sources from each field above the sig- 
nificance criterion of 10~^. Infrared (IR) counter- 
parts from the catalog of Barsony et al. (1997) 
are searched and the position is cross-correlated 
to that of the X-ray source. Systematic coordinate 
offset of the Chandra frame is then fine-tuned to 
fit the IR frame. After the offset correction, the 
relative position error (Icr) is 1'.'2. Wc find X-ray 
counterparts from 5 IR positions in a 2a error ra- 
dius (2'.'4) out of the 18 catalogued BDs and CBDs 
(GY 310, GY 31, GY 37, GY 59, and GY 326). 
The X-ray positions and relative offsets from the 
IR sources are given in Table 2. The X-ray counts 

^see http;//asc. harvard.edu/udocs/docs/swdocs/detect/html/ 



are extracted from a circle of a half radius of the 
point-spread function (PSF) around the X-ray po- 
sition. 

Since no apparent X-ray sources are found in a 
2a error radius (2'.'4) of the other 13 catalogued 
BDs and CBDs, we define a circle with a half ra- 
dius of PSF around each IR position. Then we 
manually count the X-ray photons in the circle. 
We note that a rather small source radius is se- 
lected so as to maximize the signal-to-noise (S/N) 
ratio, particularly, for faint X-ray sources. Never- 
theless, as is demonstrated in the 5 bright sources, 
the position error between IR and X-rays is always 
smaller than the source radius, because both have 
generally similar dependence on the source off-axis 
angle; both have the smaller values for the sources 
with smaller off-axis angle. Therefore most of the 
X-ray photons for the relevant BDs and CBDs, if 
any, may fall in the source circles. The X-ray flux 
thus counted are given in Table 2. 

The mean background counts are estimated 
from source-free regions of a 63-arcmin^ and a 
59-arcmin^ area in the ACIS-I fields for obs.l 
and obs.2, respectively. The soft band (0.5-2.0 
keV) background counts (in units of 10~^ counts 
arcsec"^) are 2.3 and 2.2, for obs.l and obs. 2, 
respectively, while those in the hard (2.0-9.0 keV) 
band are about 3 times larger, 6.8 and 6.5 for obs.l 
and obs. 2. The background counts in each source 
area are given in Table 2. We then separately cal- 
culate the confidence level (CL) of the X-ray de- 
tection for the soft, hard and total (0.5 -9.0 kcV) 
bands. Based on the Poisson statistics, the CL is 
defined as; 

No-l atN' 
N'=0 

where Nq and N^g are the detected and the back- 
ground counts in the source circle, respectively 
(sec e.g., Gchrcls 1986, and references therein). 
We set the detection (or the upper limit) criterion 
that the CL should be larger than 0.999 (=3.3(7) 
in any of the 3 energy band data. In Table 2, we 
show the maximum CL value for each source. For 
sources located in both fields of view (GY 141, GY 
163, and GY 202), wc also estimate the CL for the 
combined data as well as the separate data. How- 
ever, no higher CL value is obtained from all the 
sources. 
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As arc summarized in Table 2, bright X-rays 
are found from 3 CBDs (GY 31, GY 37, and GY 
59) in addition to previously reported sources, GY 
310 and GY 326 (Imanishi, Tsuboi, & Koyama 
2001). Detection of X-rays from a CBD, GY 84 
is also highly significant. Another BD, GY 141 
in obs.l, show faint X-rays with CL of >0.999 in 
the soft band. Using the logA^-logS relation by 
Mushotzky et al. (2000), we estimate the chance 
coincidence of a background source to fall in the 
source circle to be 0.6-1 %, where the uncertainty 
mainly comes from the ambiguity of the amount of 
interstellar absorption. Furthermore, background 
sources are most likely AGNs and may have harder 
spectra than the possible counterpart of GY 141. 
Wc thus conclude that the faint X-rays from this 
BD are real. 

3.2. X-ray Spectra and Luminosities 

For the brightest 4 sources — , 1 BD (GY 310) 
and 3 CBDs (GY 31, GY 59, and GY 326) — , 
we analyze the time-averaged X-ray spectra. The 
background spectra arc extracted from the same 
regions used in §3.1. We then fit the spectra with 
a thin-thermal plasma (MEKAL: Mewe, Gronen- 
schild, & van den Oord 1985) model, in which the 
metal abundance is fixed to 0.3 solar based on the 
fitting results of other Chandra X-ray sources in 
this region (Imanishi et al. 2001). The spectra and 
the best-fit parameters are shown in Figure 1 and 
Table 3. The best-fit plasma temperature (fcT) 
is in the range of 0.9-2.5 kcV. Small differences 
in the parameters of GY 310 and GY 326 from 
those in Table 1 of Imanishi et al. (2001), which 
is nevertheless unimportant, would be attributed 
to the more elaborated data-reduction process of 
this paper (see §2 in Imanishi et al. 2001). 

For the faint or non- X-ray sources, the spec- 
trum is assumed to be a thin-thermal plasma of 
2-keV temperature (the mean value for the bright 
sources) , while the absorption column is estimated 
from the empirical relation, A^h = 1.59x10^^ Ay 
cm~^ (Imanishi et al. 2001), where Ay is the vi- 
sual extinction derived from IR photometry (Wilk- 
ing et al. 1999; Gushing et al. 2000). We then 
calculate the X-ray luminosities or those of the 
99.9 % upper limits in the soft band (Table 2). 



3.3. Time Variability 

Figure 2 shows the light curves of (a) GY 310, 
(b) GY 31, (c) GY 59, and (d) GY 326 in 0.5- 
9.0 keV, where the background levels are repre- 
sented by the dashed lines. For GY 31 and GY 
59, clear flare-like events, typical to those of low- 
mass main-sequence and pre-main-sequence stars; 
a fast rise and slow decay, are seen. The flare- 
decay timescale (~10^ sec) is also the same as that 
of low-mass stars. 

GY 310 and GY 326, on the other hand, show 
no clear flare, but exhibit aperiodic variability of 
flux change by a factor of ~2 within 100 ks ex- 
posure. Since we see no large intensity variation 
in the background flux, we conclude that the vari- 
ability is highly significant. 

In order to search for spectral change during 
the flare, we separately make and fit the spectra 
of GY 31 in the flare and quiescent phases. A 
hint of spectral hardening during the flare is found, 
although it is not significant from the statistical 
point of view; kT is 2.5to.7 the first flare 

then decreases to 2.11q;5 keV in a rather quiescent 
phases. 

4. DISCUSSION 

We detect X-ray emission from 2 out of 8 BDs 
and 5 out of 10 CBDs in the p Oph cloud cores. 
The X-ray detection rate is thus 25% (BD) and 
50% (CBD). Neuhauser et al. (1999) claimed 
the X-ray detection of only 1 BD (GY202) with 
ROSAT/PSPC, in the same flelds of the present 
observations. Our present observations flnd no X- 
ray from GY202 with an upper limit of 1.2x10^* 
ergs s~^. Instead, we flnd a class I candidate 
WL 1 (No.l3, in Imanishi et al. 2001) within the 
PSPC error circle of GY202. This source emits 
X-rays with comparable luminosity to that of 
GY202 reported with ROSAT/PSPC. Conversely, 
Neuhauser et al. (1999) foimd no X-ray with an 
upper limit of ~10^^ ergs s~^ from all the Chan- 
dra detected sources. They, however, estimated 
the X-ray luminosity assuming kT = 1 keV and 
Ay = a, few mag (i.e., A^h = a few xlO^^ cm~^), 
which may be signiflcantly lower than the real case 
(see Table 3). We hence re-estimate the upper 
limit with the same assumption adopted in this pa- 
per (§3.2), then flnd that the i? 05^ T upper limit 
could be nearly one or two order magnitude higher 
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than those of the original paper. If we adopt these 
newly estimated upper limits, we find no signifi- 
cant flux change from the ROS AT to the Chandra 
era. 

In Table 2, we display the ratio of the X-ray 
{Lx) and bolometric (Lboi) luminosities. Lboi 
is roughly proportional to the area of the pho- 
tosphere, hence Lx/Lboi is a good indicator of 
coronal X-ray activity per unit surface area. The 
Lx/Lhoi value lies between 10~^-10~^, which is 
similar to low-mass pre-main-sequence stars (e.g., 
Imanishi et al. 2001) and dMe stars (Giampapa 
et al. 1996), but is significantly higher than that 
of the quiescent solar corona; '~5xl0~'' (Vaiana & 
Rosner 1978). It should also be noted that most of 
the upper limits of Lx / iboi of X-ray non-detected 
BDs and CBDs are scattered around 10~^-10~^, 
which are comparable with or slightly lower than 
those of X-ray detected BDs and CBDs. This leads 
us to suspect that these X-ray non-detected BDs 
and CBDs may emit X-rays below but near at the 
current sensitivity limit of Chandra. 

We also obtain the X-ray spectra from 1 BD and 
3 CBDs for the first time (§3.2), which are fitted 
with a thin thermal plasma model of ~2 keV tem- 
perature. Solar-like flares are also detected from 
2 CBDs (GY 31 and GY 59, §3.3). These X-ray 
features are similar to those of low-mass stars. To- 
gether with the high Lx/^boi value, we suggest 
that X-rays from BDs and CBDs are attributable 
to magnetic activities like low- mass stars. Detec- 
tion of highly variable radio emission from GY 
31, which probably arises from non-thermal gyro- 
synchrotron emission, also indicates the existence 
of strong and time variable magnetic activity (P. 
Andre, private communication; see Neuhauser et 
al. 1999). 

One debatable issue is the mechanism of mag- 
netic field amplification. Brown dwarfs are fully 
convective, hence the standard dynamo (a-u) dy- 
namo) mechanism may not work. This situation 
is the same as low-mass protostars. As one pos- 
sible scenario, Montmerle et al. (2000) proposed 
an amplification of the magnetic field between the 
star and the circumstellar disk by twisting through 
the differential rotation. If BDs and CBDs emit 
X-rays with the same mechanism. X-ray activity 
should vanish when the circumstellar disk disap- 
pears as the star evolves. Figure 3 shows the rela- 
tion between the excess emission at i^-band and 



the Lx/Lboi value. The excess emission is defined 
as r/c = Fk^^/Fk, where Fk^^ is the flux from ex- 
cess (circumstellar) emission and Fk is the stellar 
flux at K-ha,nd (Wilking et al. 1999; Gushing et 
al. 2000), hence > indicates the presence of 
the circumstellar disk. X-ray detected BDs and 
CBDs are uniformly distributed over two regions 
of Tfc > and Vk < 0; we see no hint that X-ray 
activity decreases as the circumstellar disk disap- 
pears (rfc < 0). We therefore conclude that the 
existence of the disk is not the key parameter in 
X-ray emission of BDs and CBDs. Another type 
of dynamo, such as turbulent-driven (a^) dynamo 
(Durney, De Young, & Roxburgh 1993), is con- 
ceivable. 

The authors express their thanks to Dr. Thomas 
Preibisch for critical refereeing and useful com- 
ments. The authors also acknowledge Dr. Yohko 
Tsuboi and Jun Yokogawa for useful discussions 
and comments. K.I. and M.T. are financially sup- 
ported by JSPS Research Fellowship for Young 
Scientists. 
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Fig. 1.— Spectra of (a) GY 310, (b) GY 31, (c) GY 59, and (d) GY 326. The upper panels show data points 
(crosses) and the best-fit model (solid line), while the lower panels are the data residuals from the best-fit 
model. 
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(b) GY 31 
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Fig. 2.— Light curves of (a) GY 310, (b) GY 31, (c) GY 59, and (d) GY 326 in the 0.5-9.0 keV band. 
The background level is shown by the dashed line. For GY 31, the background line (~2xl0~^ counts s~^) 
overlaps with the bottom solid line. 
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Fig. 3. — Plots of the excess emission at K-band (r^) and the luminosity ratio between X-ray and bolometric 
{Lx/L\jo\)- Circles and squares represent BDs and CBDs, respectively. X-ray detection and non-detection 
are indicated by filled and open symbols, respectively. Arrows indicate the upper limits. GY 31 and GY 163 
are not included in this figure because no rk value has been obtained. 
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Table 1 

Log of the Chandra ACIS-I Observations on the p Oph Cloud Cores 



Obs.ID 


Sequence ID 


Date 


R.A." 


Dec." 


Exposure 








(J2000) 


(J2000) 


(ks) 


1 


200060 


2000 Apr 13-14 


16''27"18!l 


-24°34'21''9 


100.6 


2 


200062 


2000 May 15-17 


16''26"35!3 


-24°23'12'.'9 


96.4 



^The position of the telescope optical axis. 
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Table 2 

BONA-FIDE AND CANDIDATE BROWN DWARFS IN THE p OpH ClOUD CORES 



Name 


Sp. type*" 


id'" R.A.'= 


Dec.^ 


offset"^ PSF" X-ray counts' 


1 - CL^ 


Lx" log(Lx/iboi) 


Ref. 






(J2000) 


(J2000) 


(") (") Soft Hard BGD 




(10^* ergs s-^) 





Brown Dwarf (BD) 



CRBR 14 


M7.5 


2 












5.7 


3 


2 


0.5 


ND (S) 




<0.55 


< 


-4.3 


1 


GY 10 


M8.5 


2 












3.2 


2 





0.2 


ND (S) 




<0.48 


< 


-4.5 


1 


GY 11 


M6.5 


2 












3.3 


1 


1 


0.2 


ND (T) 




<0.12 


< 


-3.8 


1 


CRBR 31 


M6.7 


2 












7.1 





2 


0.9 


ND (H) 




<0.56 


< 


-3.6 


2 


GY 64 


MS 


2 












3.7 


1 





0.2 


ND (S) 




<0.93 


< 


-3.5 


1 


GY 141 


MS 


1 












9.0 


7 


1 


1.5 


S.OxlO"*' 


(S) 


0.25 




-3.6 


2 






2 












22.3 


11 


23 


8.6 


ND (S) 


<0.33 


< 


-3.5 




GY 202 


M7 


1 












9.4 


2 


4 


1.6 


ND (S) 




<1.2 


< 


-3.7 


1 






2 












16.8 


3 


9 


4.9 


ND (S) 




<2.0 


< 


-3.4 




GY 310 


MS. 5 


1 


27" 


°3S!6 


38' 39'. 


'1 


1.3 


9.2 


182 


59 


1.5 


<1.0xlO"^-'' 


(T) 


9.3 




-3.4 


1 


candidate Brown Dwarf (CBD) 
































CRBR 15 


M5 


2 












4.2 


2 


1 


0.3 


ND (S) 




<1.1 


< 


-4.2 


1 


GY 5 


M7 


2 












4.7 


2 


1 


0.4 


ND (S) 




<0.16 


< 


-5.2 


1 


GY 31 


M5.5 


2 


26" 


^25^2 


23' 25' 


'7 


0.2 


2.5 


41 


350 


0.1 


<1.0xlO"^^ 


(T) 


93 




-3.5 


1 


GY 37 


M6 


2 


26" 


^27^8 


26'43'. 


'9 


0.8 


4.2 


7 


5 


0.3 


3.7x10""* 


(S) 


0.59 




-4.2 


1 


GY 59 


M6 


2 


26" 


^31=4 


25'32' 


'2 


0.7 


2.6 


17 


22 


0.1 


<1.0xlO~i^ 


(T) 


2.6 




-3.9 


1 


GY 84 


M6 


2 












1.9 


2 


2 


0.1 


1.3x10"*' 


(T) 


0.85 




-4.7 


1 


GY 163 


M2.5 


1 












12.4 


5 


12 


2.8 


ND (T) 


<12 


< 


-3.S 


1 






2 












11.9 


3 


7 


2.4 


ND (S) 




<10 


< 


-3.8 




GY 218 


M4.0 


1 












7.2 


2 


3 


0.9 


ND (S) 




<0.45 


< 


-3.6 


2 


CRBR 67 


M4.4 


1 












3.2 








0.2 


ND (-) 




<0.31 


< 


-3.5 


2 


GY 326 


M4 


1 


27° 


'42!8 


38'51'. 


'1 


1.5 


11.4 


67 


65 


2.3 


<1.0xlO"^^ 


(T) 


28 




-3.4 


1 



^Spectral type determined from the water vapor absorption. 
^Observation identification (see Table 1). 

■^Position of the identified X-ray source which is picked up with the wavdetect command (see §3.1). Right ascension and declination for all sources are at 16^ 
and -24°. 

'^Offset between the IR and the X-ray sources. 

^Radius of the point-spread function. 

*^X-ray counts in the soft (0.5—2.0 keV) and hard (2.0—9.0 keV) bands, and predicted background counts in the soft band within the source circle. Background 
counts in the hard band is ^3 times larger than those in the soft band. 

^CL is the confidence level of the source detection (see §3.1). Parentheses indicate the energy band in which the maximum CL value is obtained (S: 0.5-2.0 
keV, H: 2.0-9.0 keV, T: 0.5-9.0 keV). ND represents non X-ray detection with CL < 0.999 (i.e., 1-CL > 10"^). 

^Absorption-corrected luminosity in 0.5-9.0 keV. The distance to each source is assumed to be 145 pc (de Zeeuw et al. 1999). 



References.— (1): Wilking et al. (1999); (2): Gushing et al. (2000) 



Table 3 

Best-Fit Parameters of the Brightest 4 Sources 



Name 




log(£;M)'' 




/d.o.f 




(keV) 


(cm-=*) 


(10==cm-2) 




GY 310 


1.7(0.9-2.2) 


51.6(51.5-52.1) 


0.8(0.5-1.5) 


57.8/44 


GY 31 


2.2(1.7-2.9) 


52.6(52.4-52.8) 


5.9(5.1-7.2) 


25.1/33 


GY 59 


2.5(>1.0) 


51.0(50.6-51.7) 


1.4(0.5-3.0) 


0.556/4 


GY 326 


0.9(0.7-1.2) 


52.1(51.7-52.4) 


2.3(1.6-3.0) 


34.8/28 



^Parentheses indicate the 90% confidence limits. 
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